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We theoretically show the negative refraction existing in Mobius molecules. The negative re¬ 
fractive index is induced by the non-trivial topology of the molecules. With the Mobius boundary 
condition, the effective electromagnetic fields felt by the electron in a Mobius ring is spatially inho¬ 
mogeneous. In this regard, the Dn symmetry is broken in Mobius molecules and thus the magnetic 
response is induced through the effective magnetic field. Our hndings open up a new architecture 
for negative refractive index materials based on the non-trivial topology of Mobius molecules. 
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I. INTRODUCTION 

It is recognized that materials with simultaneously 
negative permittivity and permeability also smport the 
coherent propagation of electromagnetic fields [I| , like the 
usual case where both permittivity and permeability are 
positive. Due to a reversed phase velocity with respect 
to the group velocity 0, the boundary condition makes 
the light to bend in an unexpected direction as if the re¬ 
fractive index appears negative according to the Snell’s 
law. Hence, the effect is termed as negative refraction. 

Negative refractive index materials are very promis- 
ing^e.g., in achieving the electromagnetic field cloaking 
[1, l^, facilitating the sub-wavelength imaging and 
crime scene investigation Q. However, such materials 
do not naturally exist for the absence of magnetic re¬ 
sponse at the same frequency as electric response. To 
overcome this difficulty, delicately designed metamateri¬ 
als have been proposed and tested. Metamaterials rely 
on structural designing of the sub-wavelength unit cells 
to tune the electromagnetic resonant characteristic Q . A 
key element commonly involved in such artificial archi¬ 
tectures is the configuration of split-ring resonator [M3|- 
This structure is analogous to an LC resonator with char¬ 
acteristic frequency wq Ij^/LC with L and C being 
self-inductance and capacitance of the resonator, which 
is tunable by engineering. A metamaterial with split¬ 
ring resonators can negatively respond to the magnetic 
field when w > wo Hi- Although the negative refrac¬ 
tion has been demonstrated at the visible wavelength, 
manufacturing of microstructure at the size of 30-100nm 
is required |ll| . Such technically requesting fabrication 
makes it still an open challenge to scale the materials to 
3D bulk. 
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Instead of fabricating those delicately designed struc¬ 
tures from the conventional “top-down” approach, an at¬ 
tractive alternative is to generate the dielectric media 
with self-assembled functional atoms or molecules di¬ 
ll^. This “bottom-up” approach is particularly inter¬ 
esting due to the small size of molecules and their as¬ 
sociated quantum effects. For example, by utilizing the 
quantum interference among multi-level atoms, it is pos¬ 
sible to suppress the absorption meanwhile keep reason¬ 
able optical response We notice that the molecular 
ring can not effectively respond to the magnetic field, be¬ 
cause there does not exist an y sp lit to make the symmetry 
breaking as shown in Refs, Appendix lAl 

In this paper, instead of relying on the split-ring con¬ 
figuration, we describe a potential negative-refracting di¬ 
electric medium with the Mobius molecules [il,[il. Such 
molecules were theoretically proposed [l^ and had been 
fabricated [20l - [^ in several experiments. Compared 
with other molecular-ring-like annulenes, the symmetry 
of Mobius molecule is lowered by its boundary condition, 
i.e., the usual symmetry is broken here to C2 El- 
Similar to a superconducting ring with magnetic flux go¬ 
ing through [23j, this particular boundary condition for 
Mobius molecule could be canceled by some local unitary 
transformations. After the transformation, the motion of 
an electron in the Mobius ring is subject to an effective 
spatially-inhomogeneous magnetic field, even though the 
original system is exposed to a physically uniform elec¬ 
tromagnetic field [ 1 ^. Both electric and magnetic re¬ 
sponses can be induced at the same transitions and thus 
both negative permittivity and permeability can be ob¬ 
served at the same frequencies. In this sense, the negative 
refraction in Mobius molecules is intrinsically caused by 
the non-trivial topology of the molecules. Motivated by 
this discovery, we propose realizing a new kind of meta¬ 
materials based on the Mobius molecules. 

This paper is organized as follows. In the subsequent 
section, the Mobius molecule is briefly introduced follow¬ 
ing Ref. [ 1 ^. In Sec. Illli the allowed transitions between 
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FIG. 1: (color online) (a) Schematic of a molecular ring with 
Mobius topology: carbon atoms are shown by black circles, 
with the bonding atoms linked together, (b) Energy spectrum 
Ekcr of the molecular ring. Two energy bands are denoted by 
their different pseudo spin labels a =t and Allowed inter¬ 
band transitions from the ground state are indicated by the 
colored arrows, where dashed {x and y) and solid (x, y, and 
z) arrows indicate the possible field polarizations. 


different eigenstates of a Mobius molecule induced by ap¬ 
plied electromagnetic fields are discussed under dipole 
approximation. Then, in Sec. lIVI we investigate the possi¬ 
bility of having simultaneously negative permittivity and 
permeability in a medium consisting of non-interacting 
Mobius molecules. In Sec. El the negative refraction is 
analyzed at a planar medium interface for two different 
incident conditions. And discussions regarding loss and 
bandwidth of negative refraction are also given. Finally, 
the main results are summarized in Sec. VI. 


II. THE MOBIUS MOLECULAR RING 

The Mobius molecular ring has several physical realiza¬ 
tions, e.g., through graphene [2^ and the non-conjugated 
molecules [1^ . For a general double-ring system with 2N 
atoms, as shown in Fig. llTa), the Hiickel Hamiltonian for 
a single electron reads [l^ 1^ 


N-l 


H = 


[A]MA,-^(A]A,+i+h.c. 


( 1 ) 




where 


A, = 


M = 


e -V 
-V -e 


( 2 ) 

( 3 ) 


The difference between the Mobius molecular ring and 
the common chemical annulenes lies in the boundary con¬ 
dition [l^. Here, the Nth nucleus in the A ring is exactly 
the 0th nucleus of the B ring. Therefore, 

oo = bN, (4) 

bo = aN (5) 

imply that the operators do not obey the periodical 
boundary condition. Further calculations in Appendix 1X1 
demonstrate that the vanishing magnetic dipole for the 
perfect ring results in the absence of negative refraction. 

With a local unitary transformation Uj (see Ref. (l^ 
and Appendix IB|| . the Hamiltonian becomes 


N-l 


^ = E - C(B]QB,+i + h.c.)], (6) 

3=0 


where 


Q 


UjA, = 

■giV2 0 
0 1 


CR J ’ 


(5 = 2Tr/N, 


( 7 ) 

( 8 ) 
( 9 ) 


the Pauli matrix (Tz is written in the pseudo spin space, 
i.e., Bjo-jB^ = We emphasize that the 

new operators Cja-’s satisfy the periodical boundary con¬ 
dition. Equation represents a pseudo spin in a hc- 
titious ring [l^. For the spin-up state, the ring is fur¬ 
ther penetrated by a perpendicular homogeneous mag¬ 
netic field, as indicated by the phase factor exp{i6/2) 
in the hopping term. More insight on the non-trivial 
boundary condition stems from the corresponding energy 
spectrum, cf. Fig. [TKb) and Appendix [BJ 


c 

Ek^ = V-2^cos{k--), (10) 

Eki = -V-2^cosk, (11) 


which are obtained by performing the Fourier transform 
on Bj with the resulting transformed operator denoted 
as Cfc and k = IS {I = 0,1, - ■ ■ N— 1). Compared with the 
spectrum of a topologically trivial annulene, the effective 
magnetic flux results in a shift of the Ek'f band such that 
the two lowest states become degenerate. Similar shift¬ 
ing behavior has also been s ugg ested for a cyclic polyene 
under a real magnetic field |26j . 


Here, the fermionic operators a](&j) create a localized 
atomic-orbital \(l>j+) (|(/)j_)) at the jth nuclear site of the 
A (B) ring respectively; 2e describes the on-site energy 
difference between atoms of two rings; The inter (intra) 
-ring resonance integral is denoted by V (^). Hereafter, 
we consider a special case where all atoms are of the same 
species, i.e., e = 0 [1^, as indicated in Fig. [T](a). 


III. TRANSITION SELECTION RULES 

The general interaction Hamiltonian of a molecule with 
external electromagnetic field is complicated. But due to 
the small size of a molecule (typically around a few nm) 
with respect to wavelength the dipole approximation is 
applicable. In the linear response regime, non-vanishing 
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electromagnetic response implies that the corresponding 
dipole operators have off-diagonal elements between dif¬ 
ferent eigenstates of H. Hence it is useful to analyze the 
transition selection rules under the dipole approximation. 
We investigate the dipole transition selection rule of the 
Mobius molecular ring under the perturbation of an ex¬ 
ternal oscillating electric field with amplitude Eq and fre¬ 
quency uj. The interaction Hamiltonian is written under 
the dipole approximation as 


integral. For the Mobius molecular ring, /3y is either ^ 
or V, depending on whether the atoms are on the same 
ring or not. Consequently, the magnetic dipole operator 
is explicitly given as 

m = -eJ2[^](yS+-ay)A, + C(*A]N,A,+i + h.c.)] 

( 17 ) 

with the effective area 


H'^ = -d - Eq cos ujt, ( 12 ) 

where the electric dipole operator d = —ef and r is the 
position operator for the electron. 

The transition selection rule can be inferred from the 
matrix elements of H'^ between the eigenstates of El, 
which are the linear combinations of the atomic orbitals. 
In the following we ignore the overlap integrals from dif¬ 
ferent atomic orbitals and assume [l^ 

{4^js - djj/5ss' Idjs-) (^^) 

where Rj+{Rj-) denotes the position vector of the jth 
nuclear site at the A(B) ring. For the Mobius molecular 
ring with radius R and width AW, the nuclear positions 
are explicitly written in the molecular coordinate system, 
cf. Fig. [TKa) and Ref. [l^, as 

Rj± = (i? ± VF sin -^) cos (pjix -|- (i? ± VF sin —) sin (fijCy 
±W cos^iz, ( 14 ) 


where 


(Pj = j6 ( 15 ) 

is a rotating angle which locates atoms on the rings. 

The transition selection rules for the electric dipole op¬ 
erator are briefly summarized as, cf. Appendix ICl 

|fc, t), \k,i) "4 \k + 2(5, t), 

\k,a) |fc ± (5, cr'), |fc, 4.) ^ + <5, t) ■ 

Here, \k, a) denotes an eigenstate of H with energy Eka- 
The superscripts x, y, z indicate the electric field polariza¬ 
tions. Those allowed inter-band transitions are schemat¬ 
ically shown in Fig.[I](b). 

Since the negative refraction depends on the simulta¬ 
neously negative permittivity and permeability, we fur¬ 
ther investigate the transition selection rules for the mag¬ 
netic dipole rfi. For molecular systems the current is re¬ 
stricted to the chemical bond, thus rh is written as 
TO = JijSij, where Sij is an effective area specified 

below. The bond current Rj flows from the ith atom to 
the jth atom, and by the tight-binding approximation 
the summation only runs over bonded-atoms pairs. The 
operator corresponding to the bond current is given by 
Jij = ie/dijolaj + h.c. [l^, where Rj is the resonance 


Stf ^ \r,^ X Ryp 

{a, (3 = ±) and the vectorial matrix 





0 

1 


( 18 ) 


( 19 ) 


Under the dipole approximation, the interaction 
Hamiltonian describing the coupling of the Mobius 
molecular ring to an external oscillating magnetic field 
with amplitude Bq and frequency oj is similar to Eq. m- 


H'g = —TO • Bq cos wt. 


( 20 ) 


The explicit expressions of the matrix elements are com¬ 
plicated in the case of magnetic perturbation and thus we 
list them in Appendix [CJ Through straightforward cal¬ 
culations, we find the following transition selection rules 
as illustrated in Fig.[T](b): 

xvz xy 

|fc,;)v^|fc,t), |fc,;)v^|A: + 2(5,t), 

xy xyz 

l^i t) ^ 4-) j 1^14-) 1^ + t) I 

where only the inter-band transitions are shown since in 
the limit of large N, the spectra EkRs become continu¬ 
ous with respect to k and thus the inter-band transitions 
are more relevant to the response at high frequency as 
compared to the intra-band transitions. A comparison 
between Eq. m and Eq. (HU) indicates that the selec¬ 
tion rules for the inter-band transitions are almost the 
same for both electric and magnetic couplings. 


IV. PERMITTIVITY AND PERMEABILITY 

In order to realize the negative refraction, the simul¬ 
taneously negative permittivity ^ and permeability ^ 
are required [l[. To evaluate those quantities, let us con¬ 
sider a negative index medium realized by single crystal 
of Mobius molecules, with the same orientation for each 
molecule. Although the permittivity and permeability 
are originally treated as scalars, for anisotropic medium 
they could generally be second order tensors that depend 
on the frequency w of external fields E and H. ^ and ^ 
are calculated by considering that the electric displace¬ 
ment D and magnetic induction B are related to the 
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polarization P and the magnetization M of the medium 


D — E — SqE P, (^2) 

B = fio'^ H = noH + fioM. (23) 


Here, P and M are quantum mechanical averaging of 
the electric and magnetic dipole operators in the per¬ 
turbed molecular ground state. In the linear response 
regime, applying the Green-Kubo formula yields [30j, cf. 
Appendix [Pl 


——Re l ^g,ka{dkrj,g ' E) 

hvo ^ 10 - Aka + *7 

k^<7 

Mo f '^g.,k<7 ijk^ka^g ' 

OJ — Aka + ij 

fe,cr 


(24) 

(25) 


Here, Aka = Eka — Eg is the molecular resonant transi¬ 
tion frequency between the excited and ground states of 
Hamiltonian H. Og^ka = {9\0\k,a) is the matrix element 
of operator O between the unperturbed molecular ground 
state \g) and excited state \k,a). The ground state is 
excluded in the summation, as indicated by the prime. 
The lifetimes r = I /7 of the excited states are assumed 
to be identical, vo is the volnme occnpied by a Mobius 
molecule in the medium. In the above calculation, the 
inter-molecule interactions have been neglected, i.e., the 
polarization and magnetization are assumed to be the 
contribution by a collection of non-interacting particles. 
This approximation is valid either in the low density limit 
or for medium with crystal structure, e.g., Zinc-Blende 
[ 2 ^. Also, we have verified that both the central fre¬ 
quency and bandwidth of negative refraction are not sub¬ 
stantially modified when the inter-molecular interaction 
is considered at the level of Lorentz local field theory [ 2 ^ 
as proven in Appendix [El 

Because ~ 1 — 0(|P|/|if|) and — 

0(|M|/|i?|), Eqs. (I24I25I) imply that the negative val¬ 
ues of ^ and ^ occur necessarily around the molec¬ 
ular resonant transition frequencies Aka- Furthermore, 
the responses exist only if the corresponding transitions 
are allowed by both electric and magnetic dipole inter¬ 
actions with fields. As for the Mobius molecular ring, 
it is particularly interesting that the inter-band transi¬ 
tion |fc,t) \k',i) could be allowed both electrically 
and magnetically by the dipole couplings. Besides, the 
typical value for resonance integral is around the order 
of a few eV [HI, , already in the region of visible fre¬ 
quency. Therefore, simultaneously negative permittivity 
and permeability in the visible frequency regime might 
be realized in a medium containing molecules with the 
Mobius topology. 

To quantitatively investigate this possibility, ^ is ex¬ 
plicitly calculated for the Mobius ring by using the matrix 
elements dka.g- For positive V and the ground state 
of H is 1 0 , 4 .) and thus Eg = When co is near the 

inter-band transition frequency Aof and for sufficiently 




FIG. 2: (color online) Relative permittivity and perme¬ 
ability jit as a function of the detuning Acu = lo — Ao,.|. 
between the frequencies of incident field to and transition 
|0i4) ^ |0it)- Only negative diagonal elements of ^ and 
in their corresponding principal-axis coordinate systems are 
shown. Here we adopt following parameters V = f = 3.6eV 
[3l| . W = 0.077nm [^, R — NW/tv, 7 “^ = 4ns [3^ . 


long lifetime r, ^ is simplified from Fq. (1^^ by only 
including the summation term whose denominator con¬ 
tains a transition frequency that is equal to Ag-t-. With 
this approximation, the relative permittivity is simplified 
in the molecular coordinate system as 




1 - 77 '(w) 0 0 

0 i-VM -27'(w) 

0 —2ifj'{uj) 1 — 4?7'(a;) 


(26) 


where 


Tj{u}) 


1 _ e^IF^ _ 

Sheovo Lo — 2V — 2^(1 — cos f) + * 7 ' 


(27) 


and rj' {rj") is real (imaginary) part of rj. The permittivity 
tensor is not diagonal in the molecular coordinate system 
because of the non-zero off-diagonal term Since 

the tensor is real-symmetric, ^ is diagonalized in the 
principal-axis coordinate system by a rotation around the 
X axis. This yields 


£1 = 1- 577 ', 


(28) 


1 — 77 ' and 1 for relative permittivity along three principal 
axes, respectively. In Fig. [5Ka), the relative permittivity 
is numerically demonstrated with the following parame¬ 
ters, V = ^ = 3.6eV [m, IT = 0.077nm (H, ttR = AW, 
7 “^ = 4ns [ 3 ^, which are obtained by fitting the spectra 
in experiments. 

On the other hand, the tensor of the relative perme¬ 
ability ^ is represented by 


)lt(ui) = 


1 — a^rj' 0 0 

0 1 — a ^ rj ' —20/377' , 

0 —20/377' 1 — 4/3^77' 


where 


a 

p 


^[T-b^(cos,5-cos|)], 




(29) 


(30) 

(31) 













5 


The permeability tensor is not diagonal in the molecular 
coordinate system either. Because a and /3 are generally 
different from 1, and ^ cannot be diagonalized by 
the same rotating transformation. In other words, the 
principal axes in which is diagonal do not coincide 
with the principal axes of ^. This yields 

m = !-(«'+4^')^', (32) 

1 — and 1 for relative permeability along three prin¬ 
cipal axes, cf. Fig. [5Db), respectively. 


V. NEGATIVE REERACTION AT MEDIUM 
INTERFACE 




FIG. 3: (color online) Schematic plots of the reflection and the 
refraction at the medium interface z = 0. (a) The electric field 
and (b) the magnetic field of the refracted light is polarized 
along the ^-direction, respectively. 


The Mobius medium is anisotropic in the sense that the 
relative permittivity ^ and permeability are tensors 
rather than scalars. Therefore, apart from demonstrat¬ 
ing simultaneously < 0 and ^ < 0 in some frequency 
region, the more concrete way to show the existence of 
negative refraction is to investigate the behavior of re¬ 
fracted electromagnetic waves at the medium interface. 
Here, we investigate the reflection and refraction at a 
planar interface between the Mobius medium and the air 
for two specific incident configurations. The results show 
that the medium is “left-handed” for E-polarized propa¬ 
gating mode [ 3 ^. 

The behavior of a propagating wave inside a medium 
is captured by its phase and group velocities. For ex¬ 
ample, when the negative refraction was first introduced 
in Veselago’s seminal paper [I|, a “left-handed” material 
was described as a medium in which the electromagnetic 
wave propagates with the opposite phase velocity with 
respect to the group velocity. In media where the permit¬ 
tivity and permeability tensors are symmetric, the group 
velo city is along the same direction as the Poynting vec¬ 
tor [^. In such case, the reversal of the phase velocity 
to the Poynting vector has been applied as a criterion for 
“left-handed” materials: kt ■ St < 0 [s^, where kt and 
St are respectively the wave vector and Poynting vector 
of refracted field. Besides this, if we consider the refrac¬ 
tion from the medium interface, the causality requires 
that the normal component of its Poynting vector is in 
the same direction as that for the incident light [s^. In 
summary, we adopt the following three criteria as charac¬ 
teristics of a “left-handed” medium: (i) The wave vector 
kt in the medium is a real-valued vector; (ii) The wave 
vector kt and Poynting vector St in the medium satisfy 
kt ■ St < 0; (iii) The normal component of the Poynting 
vector remains the same sign across the medium inter¬ 
face, e.g., if interface normal is Cz then SizStz > 0, where 
Si is the Poynting vector of the incident light. 

Let us consider a linearly polarized monochromatic 
light which is incident from the air onto the interface 
of Mobius medium as illustrated in Fig. |31 Let z = 0 be 
the medium interface and suppose that the incident plane 
is the y-z plane and 9 is the angle between the incident 


wave vector and the z axis. We consider two independent 
incident configurations respectively: 

Ei = {E-polarized) (33) 

and 

Hi = {H-polarized) (34) 

where ki = kiyiy kizCz is the wave vector of incident 
light. These two configurations are known as E-polarized 
and H-polarized respectively, as the electric and magnetic 
fields of the refracted light are perpendicular to the re¬ 
fracted wave vector [^, respectively. 

A. E-polarized Incident Configuration 

For the specific incident configuration given by 

Eq. (1551) . the three criteria for a “left-handed” medium 
could be checked one by one. It follows from the bound- 
ary conditions derived from the Maxwell’s equations that 

Cz X {Ei -\- Ej. — Et) = 0, 
ez X {H, + Hr- Ht) = 0, (35) 

kiy — ^ty — ki sin^, 

where the subscript r denotes the reflected field back to 
the air. In accordance with above boundary conditions, 
the electric fields of the reflected and refracted lights 
could be written as 

Er = (36) 

Et = (37) 

where we also use r and t to represent the reflection and 
refraction coefficients respectively, and 1 r = t. To 
determine the refracted wave vector, we combine the two 
Maxwell’s equations kt Et = ^fJ-ot^Ht and kt ^ Ht = 
—Lueo^Et- by multiplying both hand sides of the first 
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equation with ^ ^ and then inserting it to the second 
equation, the following equation could be derived, i.e. 

kt X [{^lt)~^{kt X Et)] = -co^HoSo^Et. (38) 

Since the equation is homogeneous in Etx, Ety and Etz, 
the necessary condition for non-vanishing refracted light 
is that the determinant of its coefficient matrix is zero. 
Combined with Et = EtxSx and ktx = 0, this necessary 
condition could be explicitly written as 


+ 2^^<f^h.yktz + = 0. 

(39) 

Inserting Eqs. (I26I29I) into the above equation yields 

ktz = {±ki^J^j.1 [(1 - T]') (1 - 4 / 32 ?y') - sin^e] 

-ySr^k,y}lix^r^. (40) 


The time-averaged Poynting vector of refracted light 
St = Re{Et X i?t*}/2 reads 

St = StyCy + Stz&Z, (41) 

where 

According to Eqs. (j29l40l4ip . for the E-polarized con¬ 
figuration, the three criteria for “left-handed” medium 
are specified as kt ■ St < 0, and 

(1 — 77 ') (1 — 4 / 3 ^ 77 ') > sin^ 6, if / 7 i > 0, (44) 

(1 — 77 ') (1 — 4^13^1]') < sin^ 0, if fii < 0, (45) 

and 

Stz = -7^^ — [-2al3r]'k^y -I- (1 - 4:[3'^r]')ktz] > 0 . (46) 

FigurellDa) shows the “phase diagram” in the O-oj plane 
as determined by the above three criteria. Obviously, 
there is a frequency window in which the applied electro¬ 
magnetic fields can be negatively refracted, while in most 
regions of the 6-uj plane the refracted light are “right- 
handed”. Besides these, there is also a considerable re¬ 
gion where the light will be totally reflected. The inset of 
Fig-HKa) also shows fii near the inter-band transition fre¬ 
quency of the individual Mdbius molecule. A comparison 
shows that there is a correspondence between the “left- 
handed” phase boundary and the zeros of 771 , which will 
be illustrated later. 

To quantitatively understand the “phase diagram”, we 
rely on the geometry relation between the Poynting vec¬ 
tor and the wave vector surface [s^. The wave vec¬ 
tor surface is defined by solutions to Eq. (1551) in the 


{ktx,kty,ktz) space (or equivalently in the {ntx,nty,ntz) 
space by fit = kt/uj-^/JIoeo). This geometry relation then 
asserts: in the medium with symmetric and the 
Poynting vector for a given k is either parallel or anti¬ 
parallel to the normal vector of the wave vector surface, 
as proven in Appendix [El 

The shape of the wave vector surface, as shown in 
Fig. [SJ strongly depends on the incident frequency. And 
those different shapes will in turn give rise to contrast¬ 
ing propagation properties of the refracted fields. Let us 
first consider an incident field with frequency lu <C Ao-f-. 
It follows from Eq. (EH) that in this case 77 ' < 0 and | 77 '| 
is of the order of unity, while a 1 and /3 <C 1. There¬ 
fore, after ignoring the second order terms of a and j3 in 
Eq. ((551) . the wave vector surface is nearly of circle shape, 
i.e.. 


~ l-’l'- (47) 

Obviously, a real solution of ktz is admitted for all inci¬ 
dent angles, cf. Fig. [5] Since <C 1 and meanwhile 
Pi , , and pr^^'^ are of order one when u <C Aq-i-, it fol¬ 

lows from Eqs. (I41M3I) that the Poynting vector St is ap¬ 
proximately along the direction of kt . Because kf St < 0 
is violated, the refracted field is “right-handed”. 

When the frequency of the incident field is increased 
to just above the lowest inter-band transition frequency 
w > Aq-i-, since rj' ^ 1, pi in Eq. (l32l) could be negative. 
By linearly combining nty and ritz to eliminate the cross 
term in Eq. (1551) , the wave vector surface is a hyperboloid, 
i.e.. 




Hy 


(1 - r]')pi 77 ' - 1 


= I, 


(48) 


where ntz=ntz sin (p + nty cos <(), nty=ntz cos p — nty sin p, 
and the mixing angle p = tan“^[—4a/3/(/3^ — a^)]/2. 
Similar non-closed wave vector surface has been sug¬ 
gested in uniaxial left-handed materials [s^. As pi < 0, 
Eq. ((4ip implies that the Poynting vector is opposite to 
the normal vector of the wave vector surface and thus 
kt ■ St < 0 , showing that the wave propagating in the 
Mobius medium is now “left-handed”. Furthermore, in 
order to fulfill the requirement for causality, the refracted 
wave vector should be on the lower branch, cf. Fig. [5] 
Equation (H51) also indicates the frequency region 
where “left-handed” propagating wave is allowed in the 
Mobius medium. Since a and /3 are small quantities, 
the sign of pi will change along with the increase of 77 b 
This sign change will make the shape of the wave vector 
surface deform from a hyperboloid to a spheroid. How¬ 
ever, because 1 — 77 ' is still negative and 771 > 0, Eq. (1551) 
will not have solution for real ntz ■ An imaginary ntz indi¬ 
cates that the wave can not propagate inside the medium 
and the incident wave is totally reflected. Therefore, the 
frequency region for the propagating “left-handed” E- 
polarized wave in the Mobius medium is determined by 
the sign change of 771 , cf. Fig. SDa), and the resulting 
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FIG. 4: (color online) “Phase diagram” of left-handedness in the 6-u) plane for (a) E-polarized incident field, (b) H-polarized 
incident field. Red and yellow regions indicate respectively the parameter regions where the refracted wave is “right-handed” 
(RH) and the incident field is totally reflected (TR) back to the air. The refracted field is “left-handed” (LH) in the blue region. 
The inserts show a magnified part in the main plot where oj is close to the lowest inter-band transition. The red dashed curve 
shows |/ii| as a function of uj near the inter-band transition. Here we use the same parameters as in Fig. [2] 



FIG. 5: (color online) Cross sections of the wave vector sur¬ 
faces of the refracted fields at different incident frequencies: 
uj <C Ao-|- (thick red, magnified by 100), ui > Aot- (thick 
blue), UJ ^ Ao-i- (thin red). The green circle, magnified by 
100, depicts the wave vector surface of the incident field {E- 
polarized). In each case, the refracted wave vector kt and the 
direction of corresponding Poynting vector St (black arrow) 
are explicitly shown, while the wave vector of the incident 
field ki = ijjy/p,QUjQ{sm6ey -\- cosOiz) is depicted as the green 
arrow. The black dashed line represents the constraint from 
the boundary condition. Here we use the same parameters as 
in Fig. [2] 


bandwidth, i.e., width of the frequency region support¬ 
ing negative refraction, is 13 = |a;i — a; 2 |, where tOjS are 
the two solutions to /ri(w_,) = 0. 

If we further tune the frequency w to be far bigger than 
Aoi, T]' approaches zero, i.e. r]' — )> 0+. Equation (HTl) 
is again valid in this case, but the radius of the wave 
vector surface for the refracted field is smaller than that 
of the incident field. However, kt ■ St > 0 suggests that 
the refracted light is also “right-handed” in this case. 


Furthermore, the total reflection of the incident light will 
happened as there is no real solution for ritz at large 
incident angle. 


B. H-polarized Incident Configuration 


The analysis for the H-polarized case is similar to that 
for the E-polarized case. In accordance with Eqs. (I35I34L 
the reflected and refracted magnetic helds are written as, 
cf. Fig. [3Kb) 

(49) 

Ht = (50) 


The equation for determining the refracted wave vector 
component ktz is similar to Eq. (1551) . i.e., 


kt X 


(tt) ^ (kt X Ht^ = -uj’^pLoEo^Ht. (51) 


By solving this equation, we obtain the refracted wave 
vector as 


ktz = {±h\Jei[{l - ri'a'^){l - dr/') - sin^ 6] 

(52) 


The time-averaged Poynting vector for the refracted 
light is St = StyCy -I- Stz&z, with 

Sty = + hye^r^), (53) 

ZLOEqSi 

and 

Stz = + ktzsi^^^). (54) 

ZCJSoSl 

The three criteria for “left-handed” medium are then 
summarized as kt ■ St <0, and 

(l — ry'a^) (1 — dry') > sin^ 0, if 5ry' < 1, (55) 

(l — ry'o:^) (1 — dry') < sin^ 0, if 5ry'> 1, (56) 




















and 


(1 - 57j')[-27^'hy + (1 - 47?')**.] > 0. (57) 

The “phase diagram” in the 9-uj plane for the H- 
polarized case is deduced according to above three cri¬ 
teria, as shown in Fig. |4Kb). Similarly to the E-polarized 
case, there are regions in the O-oj plane where the re¬ 
fracted propagating field is “right-handed” and regions 
where the incident field is totally reflected. However, in 
contrast to the E-polarized case, there are no frequency 
regions in the H-polarized case such that the propagating 
refracted field could be “left-handed”. 


C. Discussions 

In the previous subsections, we have discussed the neg¬ 
ative refraction regardless of loss. Generally speaking, 
there will be loss in the medium due to couplings of the 
molecules to the bath. The condition for negative re¬ 
fraction in the presence of loss is slightly different from 
that without loss [13, Because of the loss effect, 
the imaginary part rj" should also be considered in the 
relative permittivity and permeability. This results in 
the replacement of p' hy ij = rj' if]" in Eqs. (j26l29l) . 
Secondly, the refracted wave vector as determined from 
the Maxwell’s equations is now generally complex, and 
the phase velocit y is along the real part of the refracted 
wave vector [33, [l^. Through considering the normal 
incidence of an E-polarized field onto the medium inter¬ 
face, we can prove that the frequency region of negative 
refraction is not qualitatively changed when the loss is 
taken into consideration. 

On the other hand, as depicted in Eq. dill), the nega¬ 
tive refraction will apparently disappear if the decay rate 
of excited stats is sufficiently large. A limitation on the 
lifetime of the excited state for achieving negative refrac¬ 
tion could be deduced from the bandwidth B = |wi — a; 2 | 
according to zeros of pLi{uj), i.e. 


B = Re\ 


e2 (a2 -F 4/32) 


1 2 


Seovoh 


— 472. 


(58) 


This indicates a restriction on the excited state lifetime 

16 soVoh 


e2 (a2 + 4/32) 1^2 ’ 


(59) 


above which the negative refraction from the Mobius 
medium is expectable. Numerical verification has been 
performed and shows that the bandwidth is not qual¬ 
itatively changed with the inclusion of the loss effect. 
In experiments, it is possible to synthesize a Mobius 
ring of carbon atoms with N = 12 and radius 0.29nm 
^ , [ 22 1. And it was theoretically predicted that Mobius 
molecules with more than 60 atoms are as stable as 
Hiickel molecules [4lj. By taking the value 3.6eV for 
V and ^ [m, one finds that an excited state lifetime of 


0.51ns is enough for observing negative refraction. Since 
the excited-state lifetime of Mobius systems can reach 
as long as 350ns [H, [d^, it is reasonable to expect a 
medium with Mobius molecules as a potential material 
to show the negative refraction. 


Although in our calculation only one electron is con¬ 
sidered, our result is consistent with the more realistic 
case when all tt electrons from all atoms in the Mobius 
molecule are taken into account. In that case, when the 
spin degree of freedom is considered, two electrons with 
different spin states can stay in the same energy eigen¬ 
state |fc, cr). For the ground state of the total system 
including all electrons, all the states of the lower energy 
band will be filled. Theoretically, there could be nega¬ 
tive refraction around 4iV possible transition frequencies 
if the excited-state lifetime is sufhciently long, as implied 
by Eqs. (116121158^ . The inter-band transition in our cal¬ 
culation is just the special case of 4A^ possible transitions. 
In this sense, our simplified calculation clearly illustrates 
the key factors for demonstrating negative refraction in 
Mobius molecules. 


Eurthermore, the Hamiltonian describing a molecule 
interacting with electromagnetic field is approximated as 
dipole interaction in our calculation, cf. Eqs. (I12I20I) . 
Generally speaking, there are multi-pole contributions to 
Goulomb interaction between molecule and electromag¬ 
netic field. As long as the molecule is small, the dipole 
approximation is valid and has been freq uently used in 
the investigation of metamaterials [3,13 ■ On the other 
hand, for the sake of simplicity, the inter-molecular inter¬ 
action has been neglected in our calculation. As shown in 
the Appendix |E1 based on the Lorentz local field theory, 
both the central frequency and bandwidth of negative re¬ 
fraction have not been substantially modified when the 
interaction between molecules is taken into account. As a 
result, by modeling Mobius molecules as non-interacting 
particles, the key factors influencing negative refraction 
can be captured. 


The Mobius molecule discussed in this paper is termed 
an equilateral Mobius strip as the twist density is 
the same everywhere. In Ref. 3, it was reported 
that two conformations of tetrahydroxymethylethylene 
Mobius molecule with chemical formula C 42 H 72 O 18 have 
been synthesized. Although the equilateral Mobius 
molecule has not yet been synthesized, it was predicted 
that it could be achievable [3- On th® other hand, we 
notice that the Hiickel model with empirical parameters 
has been successfully applied to describing experimental 
data of more than 60 organic molecules with maximum 
deviation no more than 15% [U [H, . Therefore, it is 

reasonable to expect that the theoretical predictions for 
Mobius molecules could also be observed in the future 
experiments. 
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VI. CONCLUSION 

We have explored the Mobius molecular ring as a po¬ 
tential candidate for negative refraction. The previous 
investigations with the functional atoms or molecules rely 
on the conceptual analogy of the split-ring resonator for 
magnetic response [l^ . while for Mobius ring this is in¬ 
duced by its non-trivial boundary condition. Our results 
demonstrate that engineering on the topology is benefi¬ 
cial in realizing the high frequency magnetic response at 
the molecular level. This finding opens up an alternative 
approach to design molecular negative index materials, 
which is promising in achieving 3D bulk negative refrac¬ 
tion at the visible wavelength. 

We further remark that our proposal is complementary 
to the previous experimental investigation [lOj, where 
the classical metamaterial was fabricated with Mobius 
topology in a “top-down” fashion. In order to induce 
the magnetic response, their element is still based on the 
configuration of the split-ring resonator. Moreover, due 
to quantum effect, our architecture is two order smaller 
in size than theirs. 

This work was supported by the National Natu¬ 
ral Science Foundation of China (Grant No. 11121403 
and No. 11505007), the National 973-program (Grant 
No. 2012CB922104 and No. 2014CB921403), and the 
Youth Scholars Program of Beijing Normal Univer¬ 
sity (Grant No. 2014NT28), and the Open Research 
Fund Program of the State Key Laboratory of Low- 
Dimensional Quantum Physics, Tsinghua University 
Grant No. KF201502. 


Appendix A: Absence of Magnetic Dipole Transition 
in a Perfect Ring 


Gonsider a molecular ring formed by N identical atoms 
with the nearest neighbour hopping strength ^ and site 
energy e. The single electron Hamiltonian is written as 


H = 


N-l 


j=0 



(Al) 


where periodical boundary condition ag = is as¬ 
sumed. The operator creates an excitation at the jth 

atom, which is located at Rj = RcosipjCx + Rsimpjey 
and (fj is defined in Eq. (ITSl) . R denotes the radius of the 
molecular ring. 

Based on the bond current formalism in Ref. [l^, the 
magnetic dipole operator rh reads 


mx = my=0, rriz = —e^R.'^ +h.c. (A2) 

3 

Equations (|A1IA2I) indicate that a perfect ring does not 
couple to the magnetic field m. Because [m, H] = 0 
and the interaction Hamiltonian is proportional to to. 


the interaction does not mix different eigenstates of H. 
Gonsequently, the negative refraction is absent in this 
situation. 

Eurthermore, we can show that for a common double 
ring, i.e., a chemical annulene with the periodical bound¬ 
ary condition ag = a n and bg = the magnetic-dipole 
transition is not at the same frequency as the electric- 
dipole transition. As a result, there will not be the neg¬ 
ative refraction either. 


Appendix B: Diagonalization of Mobius Hamiltonian 


Before evaluating matrix elements for dipole operators, 
we solve explicitly the energy spectrum and the molecular 
eigenstates. Consider a Mobius molecular ring with 2N 
sites, which is described by the Hiickel Hamiltonian H as 
illustrated in Eq. ©■ The Hamiltonian H is in diagonal 
form when expressed in terms of and C^, i.e., 




(Bl) 


with 


Cfc = 


Et = 


dk^ 

dk]. 


N-l 


j=0 


3 > 


V — 2^ cos(fc — |) 
0 


0 

—V — 2^cosA: 


(B2) 


■ (B3) 


Therefore, the single-electron molecular eigenstates jfc, a) 
are 


|fc,t)=4tl0)- Ifc-^)=4j0). (B4) 

Here, |0) denotes the vacuum state and k = 16 with I = 
0,1,...,A-1 and5 = 27r/V. 

It is useful to express |fc, cr) in terms of the localized 
atomic-orbitals for later use. To achieve this one notice 
that, according to Ref. [l^, Bj is related to Aj by a local 
unitary transformation, i.e. 


B,= 


Cj+ 

Ci- 


= ^3^3 


(B5) 


with 


V2 


1 1 


(B6) 


and j = 0,1,..., V — 1. It follows from Eqs. (IB2IB5I) that 
the eigenstates \k,<7) can be written in terms of Aj = 
[ aj bj as 


\k,i) 


N-l 


= E 

3=0 


o-ikj 

y/N 


cU 




| 0 ) 




= E 

3=0 


y2fV 


’ g»¥>j/2 


[41 

1 

1 



|0).(B7) 
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Furthermore, by noticing that the atomic orbitals \4>j±) 
are created by acting a] or 6^ on the vacuum state, the 
molecular eigenstates are rewritten as 

\kA) = (B8) 

j=o 

N-l 

|A:,;) = + 1 '^^-)) • (B9) 

3=0 

Especially, if the resonance integrals V and ^ are positive, 
the molecular ground state is 


where d = — er is the electric dipole operator with — e 
being the electric charge and r being the position vector 
of electron, a± = {ox ± i<Ty)l2 with Ua (a = x,y,z) 
being Pauli operators in the pseudo spin space spanned 
by \k,a) and \k',a), e.g., for Eq. (ICll) is defined as 


| 5 ) = | 04 ) 


\/m 


N-l 

E 

i=o 


i\4>j+) + \(t)j-)) ■ 


(BIO) 


a, = \kA){kA\-\k,i){k,i\, (C4) 


Appendix C: Matrix Elements of Electric and 
Magnetic Dipoles 


The matrix elements of electric dipole operator are 


CfccTcI 


cjcl^s 




k±2S 


\ IV ^ „ 

Oy — {{^y 2 €z') (Xx y] (Cl) 

0^ = -e-[(ea; =F *ey) {2R + Way) 

+2VFe,a^], (C2) 

0J = -e—{Tiex - ey)azf, (C3) 


J 


ea{a = x,y,z) is the unit vector in a direction, R and 
AW are respectively the radius and width of the Mobius 
molecule. 


The matrix elements of magnetic dipole operator m = 
—ierx [H,f]/2li are summarized as follows: 


For \k,a) ^ \k,a') transitions 



e 


Ci-mC 


— |^{21F^ [cos(fc — S) — cos k] By + [lF^(cos k 
— cos{k — 26) — cos{k — 6) + cos{k + 5)) 
+4i?^(cos(fc + f) — cos(fc — |<5))]ez} 
\RW{- [V + ^(cos(fc - (5) - cos{k + f)] (cx + iBy) 
—2i^ cos I [cos(fc + jS) — cos{k — |5)] e^} 


\RW{- [V + 5(cos(fc - (5) - cos(/c + |))] {Bx - iBy) 
—2i^ cos j [cos(/c — |i5) — cos(fc + |(i)] e^} 

— ^^sinfc \W'^ sin|ey — (2i?^ + W^ cos|)sin5ez] 

(C5) 


for |/c, cr) ^ \k ±6, a') transitions 


- 0 


e 

'm 


CfcmC 


t 

k-\-6 


[cos(A; — 6) — cos(fc + (5)] {iBx + By — Bz) 

\RW{\V + ^(cos(fc + 5 ) - cos(fc - §))] {Bx - iBy) 
—*^[cos(fc — 5) — cos{k + (5) — cos(fc + |(5) 

+ cos(fc - ^)]Bz} 


— \RW{\y + ^ (cos k — cos{k + |))] ( 

—2i^ cos I [cos(fc — |5) — cos(fc + |(5)] Bz} 

\W^^ [cos(fc - I) - cos(fc + |(5)] {iBx + By - Bz) 

(C6) 
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and 


- 0 


CkmCl_g 


e 


— [cos(fc — 2 d) — cosfc] (iCx — + Cz) 


\RW{\y + ^(cosfc - cos(fc - |d))] {Cx + icy) 
+i^[cos(/c — |d) + cos(fc — 2 d) — cosfc 
-cos(fc + |)]ez} 


— jRW {V + ^ [cos(fc — d) — cos(fc — I)]} {cx + icy) [cos(fc — |d) — cos(fc + |)] {—icx + — e^) 


(C7) 


for |fc,fT) ^ |fc± 2 d, ct') transitions 


-(0 
e 


it 

^k+2S 


[cos fc — cos(fc + d)] [ex — iCy) 0 

jRW {y + ^ [cos(fc + d) — cos(fc + |)] } [ex — iey) [cos(fc + f) — cos(fc + |d)] [ex — iiy) 


(C 8 ) 


and 


- (0 

e 


CkmCl_2g 


— [cos(fc — 2d) — cos(fc — d)] [ex + iey) jRW {V + ^ [cos(fc — d) — cos(fc — |d)] } [ex + iey) 
0 —[cos(fc — |d) — cos(fc — |)] [e-x + iey) 


(C9) 


where V and ^ are respectively the inter-ring and intra¬ 
ring resonance integrals, fc is the momentum of the eigen¬ 
state, d = 2tt/N . 


Appendix D: Polarization and Magnetization 

In order to calculate the polarization of medium, we 
first consider the electric dipole [d) of a single molecule. 
According to the linear response theory 

{d) = J (Dl) 

where E[uji) is the Fourier transform of the electric field 

/ OO 

dtE[t)E^^\ (D2) 

-OO 

and the linear response function in the frequency domain 
is 

5'(i)(a;i) =-J(wi) - J*(-wi). (D3) 

Here, J(a;i) is the dipole-dipole correlation function [d^ 

poo 

J[uji) = -i dtTr[d[t)dpo]E‘^^* (Dd) 
^0 

with po = \g)[g\ and 

d[t) = (D5) 

For a monochromatic continuous driving with fre¬ 

quency uj and time-independent envelope Eq, the electric 


I- 

field is 

E[t) = Eq cos ujt, (D 6 ) 

and [d) is rewritten as 

[d) = J ^5(i)(a;i)7r.Bo[d(a;i-a;) + d(a;i+w)]e-*‘"i* 

= + S^^\-uj)e^yEo 

= + J*(-c^)]e-“‘ + [J(-a;) + J*(a.)]e“‘} 


xEq. 

According to Eq. (iDdl) . we have 

poo 

J[u}) = -i dt[g\d[t)d\g) 


(D7) 


poo 

= -lY.' 

k,a do 


= E 


/ dg^k(xdka,g 
oj - Aka + * 7 ’ 


(D 8 ) 


where dk(j,g = [k,cr\d\g) is matrix element of the electric 
dipole operator and the decay rate 7 (related to the life¬ 
time of the excited states t by 7 = 1 /r) is introduced to 
phenomenologically describe the quantum dynamics due 
to coupling to the bath. The transition frequency from 
the ground state to excited state |fc, a) is denoted as Aka- 
The electric dipole is then written as 
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-E 

k,a 


r ^g^ka^ka^g ‘ ^0 

2ft 


1 


1 


E fdg^kadktj^g ' -ft'O { 

2 ft *' 


w - + 17 UJ + Aka + il 

1 1 




_i_i_)e*‘^*] 

W - Afco- - *7 ^ + ^ka - *7 


—e 


w - Afcg. + ^7 w - Aka - *7 


(D9) 


where we have invoked the rotating wave approximation 
and ignored the terms with u + Aka in the denominator 
(43 |. Re-expressing the above equation in terms of E{t) 
yields 


(d) = 


/ ^g,k(T^ka,g ‘ -^(0 

fe,cr 

+Im^' 


ft(a; - Afeo- + n) 

dg^ka^ka.g ' E{t 


h{ui - Aka + ij) 


(DIO) 


Since the second term on the r.h.s. of Eq. (IDlOjl is related 
to absorption and loss, we only keep the real part, i.e., 


(d) 


X ^ t dg^kadka,g ' 

^ h{uJ - Aka + * 7 ) ’ 


(Dll) 


For a medium of non-interacting molecules, because 
all molecules contribute equally to the polarization, the 
polarization is obtained from Eq. (IDllIl by multiplying 
it with the molecular number density I/vq, i.e. 


m 


1 

hvo 


Re^ 

fe,(7 


fdg_}^(jdka,g ' -^(0 
LU - Aka + il 


(D12) 


For a double-ring Mobius molecule with N carbon atoms 
in each ring, the radius is i? ~ and the width 

41+ = 4i?c, where Rc = 0.077nm [s^ is the radius of a 
single carbon atom. As a result, the volume occupied by 
a single Mobius molecule is 

Vo ~ 27r(i? + l+)2l+. (D13) 


The magnetization is derived in analogy with 
Eq. (ID12I) by noticing the similarity between and 
H'e, i-e., 


( to ) 

M(t) 


X ^ /kklg^kakklka^g ' -^(0 

^ h{uj - Aka + il) ’ 

1 \ ' ^ ,ka^k(T,g ‘ -^(0 

hvo ^ OJ - Aka + *7 

k,(7 

Mo X ^ /kkig kakklka,g ‘ dd(t) 

hvo ^ W - Aka + *7 


(D14) 


(D15) 


Here, 'mka,g = {k,a\m\g) is matrix element of the mag¬ 
netic dipole operator. 


Appendix E: Local Field Correction 


Influenced by polarization of nearby molecules, the 
actual field that is exerted on an individual molecule 
could be different from the applied external field. Con¬ 
sequently, the total field at the molecule is modified as 


Etot — E + Ei, 


(El) 


where E and E^ denote respectively the external and in¬ 
ternal fields, and the latter is due to the effect caused by 
the surrounding molecules. Although the internal field 
can be obtained by mean-field approximation as 

Erae.n = -y (E2) 


where <i)(r) is the scalar potential due to charge distribu¬ 
tion inside a volume V, the approximation fails to take 
the effect of arrangement of the closest molecules into ac¬ 
count. The volume is properly chosen to include a macro¬ 
scopic number of molecules but still such small enough 
that the dipole moments for each molecules inside are 
approximately the same. To account for the arrange¬ 
ment of nearby molecules, the internal field is evaluated 
by the microscopic contribution E^ear from all charges of 
the nearby molecules in V subtracted by the mean-field 
part. This yields 281 


(E3) 


However, for molecules arranged with higher symmetry, 
e.g., cubic lattice, or totally disordered, the contribution 
from Anear Can be neglected [28|. Furthermore, by ex¬ 
panding the scalar potential up to the dipole correction 
and assuming no net charge in the volume V, the mean- 
field part is evaluated as 


Ea 


1 


= -—E 

3eo ^ 




(E4) 


where pi denotes the induced dipole moment of Rh 
molecule inside the volume V. 

For a sufficiently weak field, the induced dipole mo¬ 
ment is proportional to the exerted field, i.e. 


Pi — ^07niolEtot: 


(E5) 
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where the molecular polarizability ymoi is deduced ac¬ 
cording to Eq. (IDllI) as 


the normal of the wave vector surface. Let n be a point 
on the wave vector surface 


(b) i ^ka) 


(E 6 ) 


As the volume V is chosen such that its enclosed piS are 
approximately the same, the polarization within the vol¬ 
ume is related to the dipole moment of a representative 
molecule pi , multiplied with the molecule density \/vq 


I 


V 



(E7) 


On the other hand, the polarization P is related to the ex¬ 
ternal field through the electric susceptibility P = eoXeE. 
By combining Eqs. (IE5IE7I) . we have 

SoXeE = —7mol(E -f ^eoXeE), (E 8 ) 
Vo 3so 

from which the electric susceptibility could be obtained. 
Furthermore, by ^ = 1-1- Xe, the relative permittivity 
with the local field correction is given by 

^ = 


where in the second line an approximation similar to the 
one for Eq. (|26)) has been used. It follows from Eq. (IE9I) 
that the corresponding eigenvalues of are given by 
(3 — 2 ry')/( 3 -|- 77 '), (3 — 10 ? 7')/(3 -I- 5??') and 1 respectively. 

Compared with the results without the local field cor¬ 
rection, the bandwidth does not change qualitatively. 
From relative permittivity, the bandwidth near the lowest 
inter-band transition frequency Ao,-!- without local field 
correction is limited by the separation between two zeros 
of 1 — 5r/'(uj) = 0, while by including the local field cor¬ 
rection, it is limited by two zeros of 1 — 10y'(oj)/3 = 0. 
Hence the bandwidth in those two cases are only differed 
by a small factor and thus the bandwidth of negative 
permittivity is not significantly modified when the local 
field effect is taken into consideration. Similar proof can 
also be applied to analysis of the bandwidth of negative 
permeability. In conclusion, the local field effect would 
not substantially modify the bandwidth of negative re¬ 
fraction. 


I + (I ~ - 5 —7mol) 


3vo 

3-2ri' Q 

3+2r)' 
3+5ri' 

0 -^ 


3-1-?7' 

0 


3-I-5J7' 


Vo 

0 

Or]' 

3+brj' 

3-7r)'^ 

3-|-5t;' 


- 7 mol 


(E9) 


Appendix F: Poynting Vector and the Wave Vector 
Surface 

In this section, we follow the treatment in Ref. to 
prove that Poynting vector is parallel or anti-parallel to 


f{nj;,ny,n^) =Q, (FI) 

and 6n be such a small change to n on the wave-vector 
surface that 5n can be viewed as a vector in the tangent 
plane at that point, i.e., 5fi- df /dn = 0. Here, dfjdn is 
the normal of the wave vector surface at the point ft. 
According to the Maxwell’s equations, we have 


k X E = (jjpotTr ■ H, (F2) 

k X H = —uj£o^ ■ E. (F3) 

In combination with k = uj^JeoPoft, these yield 

5fix E Pfix 5E = c6B, (F4) 

Sri X H + ft X SH = —cSD. (F5) 


By multiplying both sides of Eq. (IF4I1 hy H, we have 

cH-SB = H ■ {Sri X E)+H ■ {fix SE) 

= S ■ Sfi— SE ■ {fix H) 

= S ■ Sfi + cSE ■ D. (F6) 

Similarly, multiplying both sides of Eq. (lF4l) by E yields 

-cE-SD = E ■ {Sfi X H) + E ■ {fi X SH) 

= —S ■ Sfi — SH ■ {fi X E) 

= -S-Sfi-cSH-B. (F7) 

It follows from Eqs. (IF6IF7I) that 

S-Sn= ^c[H ■ SB - SH ■ B + E ■ SD - SE ■ D]. (F 8 ) 

Because ^ and are independent on ft [ 3 ^, 


H-SB = H-ti-SH 



(F9) 

For symmetric ^ and i.e. 


,,(b) _ ,.Ui) Mi) — Mi) 

H"r H'r i i 

(FIO) 


we have 

H-SB = '^H,p{f^'>SHy 
ij 

= SH-t^-H 

= SH-B. (FII) 

Similarly, we can also prove that 

E-SD = ^E.el^^^SEy 

'i-3 

= J^Ersi^'^SEj 

= SE-D. (F12) 
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Inserting Eqs. (IF11IF12I) into Fq. (jF8l) yields 
S ■ Sri = 0. 


Equation (IF13I) indicates that the Poynting vector S is 
along the normal of the wave vector surface at given point 
(F13) n. 
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